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Introduction: Cosmic collisions between terres- 
trial planets resemble somewhat the life cycle of the 
phoenix: worlds collide, are consumed in flame, and 
after the debris has cleared, shiny new worlds emerge 
aglow with possibilities. And glow they do, for they 
are molten. How brightly they glow, and for how long, 
is determined by their atmospheres and their moons. 

A reasonable initial condition on Earth after the 
Moon-forming impact is that it begins as a hot global 
magma ocean 1,2 . We therefore begin our study with the 
mantle as a liquid ocean with a surface temperature on 
the order of 3000-4000 K at a time some 100-1000 
years after the impact, by which point we can hope that 
early transients have settled down. 

A second initial condition is the presence of a sub- 
stantial atmosphere. Earth today has about 270 bars of 
H20 and 50 bars of C02 near the surface, the first 
mostly in oceans and the latter mostly in carbonate 
rocks. Similarly large reservoirs are held in the mantle. 
Thus it is reasonable to start with 100-1000 bars of 
H20 and C02 in the atmosphere. The water vapor and 
carbon dioxide would be supplemented by smaller 
amounts of CO, H2, N2, various sulfur-containing gas- 
es, and by a suite of geochemical volatiles evaporated 
from the magma 3,4 . 

A third initial condition is the existence of the 
Moon. We start the Moon with its current mass at the 
relevant Roche limit. This gives Earth its shortest first 
day and maximizes the role of tidal heating. 

A fourth initial condition is the angular momentum 
of the Earth-Moon system. In canonical models this is 
held constant, so that with the Moon at the Roche lim- 
it, the Earth starts with a 5 hour day. Recent Moon- 
forming models need loss of angular momentum from 
the Earth-Moon system 5,6 . These start Earth with a ~3 
hour day and thus imply more tidal heating. 

It is well known that an atmosphere's thermal blan- 
keting effect prevents a magma ocean from cooling 
rapidly. A string of increasingly sophisticated models 7 " 
’ ’ ’ ’ have addressed consequences of thermal blan- 
keting by steam atmospheres over Earth's cooling 
magma oceans. All of these models presume H20- 
C02 atmospheres, and most ignore the Moon. 

A different set of studies has addressed the compo- 
sition of the atmosphere over a magma ocean 14 " 16,3 ’ 4 . 
Shortly after the impact the atmosphere over the lava 
also contains volatiles such as sulfur, sodium, and 
chlorine. The geochemically enriched atmosphere is 
markedly more opaque at magma temperatures than 


H20 and C02 working alone, and thus the planet cools 
more slowly and spends more time with a lava surface 
than previous models would predict. 

For Earth we use a ruthlessly simplified model. We 
characterize the Earth by its surface temperature T s and 
an interior temperature 7). In effect 7’ is a potential 
temperature, the temperature that every parcel in an 
adiabatic mantle would have were it brought to the 
surface. We also assume that the mantle is uniform, in 
the sense that it has the same viscosity at all depths and 
that it freezes everywhere at the same time. This as- 
sumption is roughly equivalent to assuming that the 
adiabat and the melting curve are parallel (a bad as- 
sumption near the surface when the surface is cold, but 
perhaps not so bad when the mantle has not yet fully 
congealed). The uniform mantle vastly simplifies the 
discussion of topics like parameterized convection that 
can depend on a global viscosity. 

We describe silicate freezing by three tempera- 
tures: a liquidus 2)i q hotter than which the mantle is 
fully molten; a solidus r so i colder than which the man- 
tle is fully solid; and a critical temperature T crlt in be- 
tween at which the rheology of the material changes 
from that of a solid with melt percolating through it to 
a liquid carrying suspended solids 11,17 . The critical 
temperature r crit corresponds roughly to what is com- 
monly meant by melting. We take 7ji q =1800 K and 
21,01=1400 K. These are arbitrary values; our results are 
insensitive to them. For simplicity we assume that the 
melt fraction is linear in 7) between 2j; q and T so i. 

We follow previous workers by using parameter- 
ized convection to link the cooling rate to the tempera- 
ture and heat generation inside the Earth 17,1 ’. In param- 
eterized convection, heat flow is related to the size and 
viscosity of the interior and the temperature gradient 
acrossthe surface boundary layer through the Rayleigh 
number. 

Viscosity is of central importance to convection 
and tidal heating. We follow 11,17 by treating the tem- 
perature dependence of viscosity in two regimes divid- 
ed by r cm . Viscosity is a strong function of tempera- 
ture, moreover it undergoes a phase transition at T=T ait 
where it jumps by more than ten orders of magnitude. 

Cooling takes place in two regimes, the first a fast- 
cooling regime in which the geothermal heat flow is 
comparable to or bigger than the globally-averaged 
absorbed sunlight, the second a slow-cooling regime in 
which the radiative cooling is in approximate balance 
with absorbed sunlight and the heat flow is small. The 



transition between regimes, when it comes, is ab- 
rupt. 



Figure 1. Thermal evolution of Earth after a canonical 
Moon-forming impact. Results are shown for a 100 bar 
atmosphere. Internal temperature 7) (mantle), surface 
temperature T s , and effective radiating temperature T e ff 
are shown against the left hand axis. The liquidus and 
solidus are also indicated. Geothermal heat flow 
(dashed curve) is shown against the right hand axis. 
The heat flow plateaus at ~100 W/m 2 while the cooling 
Earth passes through the runaway greenhouse, a phase 
that lasts from 0.5- 1.7 Myr. The runaway greenhouse 
phase ends abruptly when the monotonically cooling 
mantle becomes too viscous to sustain the high heat 
flow, at which point the surface quickly cools to the 
equilibrium temperature under —100 bars of C02. 


Moon's orbit evolves orders of magnitude more slowly 
than in conventional models. Slow orbital evolution 
promotes capture by orbital resonances that may have 
been important in the Earth-Moon system 185 . 



Figure 2. Two parameters (red) that are thought to 
play key roles in resonance capture and lunar orbital 
evolution. The velocity da/dt that the Moon recedes 
from Earth is important to assessing resonance capture. 
For conventional fixed Q, da/dt is estimated to be 10- 
100 km/yr when the Moon encounters the evection 
resonance. At these speeds resonance capture is very 
unlikely 18 . The atmosphere reduces da/dt to 0.1 km/yr 
or less, for which resonance capture is nearly certain 18 . 
The parameter k 2 /Q figures prominently in published 
work, where it is treated as constant 5 . The distance to 
the Moon is shown for context. The particular case 
shown here assumes a canonical Moon-forming impact 
and a 100 bar atmosphere. The model does not yet 
include tidal dissipation in the Moon. 


Tidal heating by the nearby Moon is a major term 
in the early Earth's energy budget. Here we assume 
that tidal dissipation in the magma ocean can be de- 
scribed by viscosity. This is a reasonable approxima- 
tion in the congealing mantle because the mantle's vis- 
cosity would be big enough to make ordinary viscous 
dissipation a big term. Where viscous dissipation dom- 
inates, Q is a computed quantity. But at very early 
times when the mantle is very hot and effectively in- 
viscid, or at modern times when the mantlis solid, or- 
dinary viscous dissipation is predicted to be negligible 
and some other source of dissipation will likely be 
more important. 

The Moon is entwined with Earth by a negative 
feedback between thermal blanketing and tidal heating 
that comes from the temperature-dependent viscosity 
of the magma ocean. Because of this feedback, the rate 
that the Moon's orbit evolves is limited by the modest 
radiative cooling rate of Earth's atmosphere, which in 
effect tethers the Moon to the Earth. Consequently the 
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